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f~>) ■ Abstract 

We study the limit as e — > of the entropy solutions of the equation d t u e + div x [A (|, u e )] = 0. 
We prove that the sequence u e two-scale converges towards a function u(t,x,y), and u is the unique 
solution of a limit evolution problem. The remarkable point is that the limit problem is not a scalar 
conservation law, but rather a kinetic equation in which the macroscopic and microscopic variables 
are mixed. We also prove a strong convergence result in L\ oc . 

1 Introduction 

Ph 

This article is concerned with the asymptotic behavior of the sequence u e € C([0, oo), L 1 1 oc (M jv )), as the 
parameter e vanishes, where u £ is the entropy solution of the scalar conservation law 

dvF{t,x) , _J2_d_ A .^ u%t ^ =0 t >o !xe R~, (1) 



> 



dt 

i=i 



: (t = 0) = u o (^x,^y (2) 



i T'Vlp fnnrtinnQ A: = AAit n\ lit 9] f= arp aQQiimprl in Kp Y^-r»prinrlir wbprp V = TT^ 



The functions At — Ai(y,v) (y G R , uGl) are assumed to be F-periodic, where F = n^ 1 (0,T J ) 
is the unit cell, and uq is also assumed to be periodic in its second variable, 
p^j ■ Under regularity hypotheses on the flux, namely A £ W^pe^iocO^^ 1 )' an d when the initial data 

u £ (t — 0) belongs to L°°, it is known that there exists a unique entropy solution u E of the above system 
for all £ > given (see [H HH [T7l El]). The study of the homogenization of such hyperbolic scalar 
| conservation laws has been investigated by several authors, see for instance OHOlIII], and in the linear 

case |14[I15|. In dimension one, there is also an equivalence with Hamilton- Jacobi equations which allows 
to use the results of [18]. In general, the results obtained by these authors can be summarized as follows: 
there exists a function u° = u°(t, x, y) such that 

u e -u°(t,x,^ ^0 in L^ c ((0,oo) x R N ). (3) 

The function u°(t,x,y) satisfies a microscopic equation, called cell problem, and an evolution equation, 
which is a scalar conservation law in which the coefficients depend on the microscopic variable y. In 
general, there is no "decoupling" of the macroscopic variables t, x, and the microscopic variable y: the 
average of u° with respect to the variable y is not the solution of an "average" conservation law. 

To our knowledge, there are no results as soon as the dimension is strictly greater than one when 
the flux does not satisfy a structural condition of the type A(y, £) = a(y)g(£,). Here, we investigate the 
behavior of the family u e for arbitrary fluxes. We prove that ((3J) still holds in some sense which will be 
precised later on, and the function u° is a solution of a microscopic cell problem. Precisely, we prove 
that even though there is no simple evolution equation satisfied by the function u° itself, the function 

f(t,x,y,Z) = l i<u a 

is the unique solution of a linear transport equation, with a source term which is a Lagrange multiplier 
accounting for the constraints on /. This statement is reminiscent of the kinetic formulation for scalar 
conservation laws (see [19, 20, 22J, the general presentation in [23], and [8J for the heterogeneous case); this 
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is not surprising since our method of proof relies on the kinetic formulation for equation (TJ. However, in 
general, it is unclear whether u° is the solution of a scalar conservation law. Thus the kinetic formulation 
appears as the "correct" vision of the entropy solutions of (HJ, at least as far as homogenization is 
concerned. 

The rest of this introduction is devoted to the presentation of our main results. We begin with the 
description of the asymptotic problem, and then we state the convergence results in the general case. 

1.1 Description of the asymptotic evolution problem 

We first introduce the asymptotic evolution problem, for which we state an existence and uniqueness 
result; then we explain how this asymptotic problem can be understood formally. 
In the following, we set, for £ M JV+1 , 



We set a(y,£) — (ai(y,£), • • • ,aN+i(y,0) •= K Ar+1 . Notice that div y ^a(y, ^) = 0. These notations were 
introduced in 



Before giving the definition of the limit system, we recall the kinetic formulation for equation |T]), 
which was derived in [8]. Indeed, we believe it may shed some light on the limit system. Let u e be an 
entropy solution of (U). Then there exists a non-negative measure m £ £ Af 1 ((0, oo) x R JV+1 ) such that 
f e = l£ <u e( t ^ is a solution of the transport equation 



In fact, this equation was derived in [8] for the function g e (t,x,£,) — x(£, u e (t, x )), where x(£i u ) = 
lo<5<« — lu<£<o, for u, £ £ E, and under the additional assumption ajv+i(y,0) = for all y £ R N . 
However, it is easily proved, using the identity f £ — g e + l^<o, that f e satisfies Q, even when aN+i(y, 0) 
does not vanish. 

We now define the limit system, which is reminiscent of equation (J4j) : 

Definition 1. Let f £ L°°([0, oo) xR w x Y xR), u £ L°°(R N xY). We say that f is a generalized kinetic 
solution of the limit problem, with initial data lf<« , if there exists a distribution M. £ T>' per ([0,co) x 
l w x Y x R) such that f and A4 satisfy the following properties: 

1. Compact support in £: there exists a constant M > such that 



2. Microscopic equation for f: there exists a non-negative measure m £ Af 1 ((0,oo) x M. N x Y x R) 
such that f is a solution in the sense of distributions of 



aN+i(y,0 = -div y A(y,£). 



dtf + ^ d x J e + -a N+1 (^,e) d ( r = d s m £ , 

r(t = 0,x,0 = l i<uo(x ,y 



(4) 
(5) 



Supp M C [0, oo) x R N x Y x [-M, M), 
f(t,x,y,0 = l »/C<-M, 
f(t t x,y,S)=0 if£>M. 



(6) 
(7) 
(8) 



div^ (a(y, £)f(t, x, y, £)) = d 6 m, 



(9) 



and Supp m C [0, oo) x R N x Y x [-M, M]. 



3. Evolution equation: the couple (f,M) is a solution in the sense of distributions of 



N 



< d t f + J2"i(y,0d x J = M, 



(10) 



. fit = o,x,y,£) = i^ <Uo { X ,y) =■ fo(x,y,Q; 
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In other words, for any test function 4> G V per ([0, oo) x R N x Fxl). 

/ / f{t,x,y,£) < dt4>(t,x,y,£) + y~]ai(y,g)d Xi <])(t,x,y,£) > dtdxdyd^ 

JO JR N xYxS. • i 



= - (0, M) v v ,-\ l i<Uo{x , y) (j)(t = 0, x, y, £) dx dy d£. 

JR N xYxR 

4- Conditions on f: there exists a non-negative measure v e M per ([0, oo) x R N xYxM.) such that 

d 6 f = -v, (11) 
< f{t,x,y,£) < 1 almost everywhere. (12) 

And for all compact set K C M. N , 

\ I 11/00 - fo\\ LHKxYxm ds — j 0. (13) 

5. Condition on Ad: define the set 

G:={i>e Lf oc {Y x R) , 9 5 V > 0, and 3/J, G M^F x R), 3C > 0, 3a_ 6 R, 

divy j {(a , 0) = —d^fji, Supp /i C F x [— C, C], fi > 0, 
1&(y,0=a- «/C<-C}. 

T/ien /or aZZ <p G 2?([0, oo) x R ) such that <p > 0, the function Ad *t, x <p belongs to C([0, oo) x 
R N ,L 2 (Y x R)), and 

V(t,i)e[0, ro )xl N , W>g£?, [ {M* t ,x<p)(t,x,-)ip<0. (14) 

JYxR 

We now state an existence and uniqueness result for solutions of the limit problem : 
Theorem 1. Let A E W p '™ loc (Y x R). 

1. Existence: let uq G L] 0C (SL N ;C per (Y)) nL°°(R N ) such that there exists a non-negative measure 
m = mo(x,y,0 such that fo(x,V,Q = lt;<u (x,y) is a solution of 

E ^7 (°i(f.0/o) + ^ (ajv+id/.O/o) = ^ (15) 
and Supp to C R n x F x [— M,M], where M = ||u ||oo- 

Assume that there exists u\, U2 G L co {Y) such that l£< Ui is a solution of lfT5|) for i = 1,2, /or 
some non-negative measures mi, mi, and 

Ui(y)<U (x,y)<u 2 (y) for a.e. x e R N ,y e Y. (16) 

TTien there exists a generalized kinetic solution f of the limit problem (in the sense of definition 
QP, with initial data fo. 

2. "Rigidity": let u G L°°(R N x Y), and let f G L°°([0, oo) x R N x Y x R) be a generalized kinetic 
solution of the limit problem, with initial data fo = lf<«„. Then there exists a function u G 
L°°([0, oo) x R N x Y) such that 

f(t, x, y, £) = l£ <u (t,x, v ) almost everywhere. 
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3. Uniqueness and contraction principle: let Uq,Vo £ L°°(M. N x Y), and let f,g be two generalized 
kinetic solutions of the limit problem with initial data T-£< Uo and l£<„ respectively. Then there 
exists a constant C > such that for all t > 0, for all R, R' > 0, 

Wf(t)-g{t)\\mB R xYxR) < e ct+R (\\u - v \\ LHBr , xY ) + e~ R ') . (17) 

As a consequence, for all uq 6 L^CU. 1 ^ xF)(l Lj oc (EL N , C per (Y)) satisfying (jT5j) and fT6|) . there exists a 
unique generalized kinetic solution f £ L°°([0, oo) x R N x Y x R) of the limit problem. 

Remark 1. Notice that for any function v £ L°°(Y), v is an entropy solution of the cell problem 

div y A(y,v(y))=0 

if and only if there exists a non-negative measure m £ M^ er (Y x R) such that the equation 

^vA a i.V^)M<v{y)) = <%m- 

is satisfied in the sense of distributions on Y x R. 

In the case where A is divergence-free this condition becomes 

N 

^dyA a i{y,OM<v{y)) = o. 

8=1 

Indeed, in that case, v satisfies 

N 
i=l 



for some non-negative measure m such that Supp m C Y x (— M, M). Consequently, 

JL / /•« 

l-M-l 



y2d yi ( [ cn(y,w)l w<v{ y) dw) =m(y,£) > 0. 



Since the left-hand side has zero mean-value on Y for all £ £ [— M, M], we deduce that m = 0. Thus, 
in the case where the flux A is divergence free, the limit system takes a slightly simpler form: conditions 
©, CH]) become 

div y (a(y,£,)f(t,x,y,£)) = 0, 

N 

d t f + J2^(y,Od x j = M, 

i=l 

Jyxm( M *t,xV) (*, x, •) V < 0, 
W> £ i^ c (^ x R), dWyiaip) = 0, and <9 5 V > 0. 

All the other properties remain the same. 

Remark 2. Assume that the flux A is divergence-free, and set 

N 

dy. 



(18) 



Ci := {ip £ L 2 loc {Y x R), £ A (ai(y, 0^(y, 0) = 0}, 



C 2 := W £ x R), drf > 0}. 

T/ien Ci, C*2 are convex sets of the vector space L^ oc (Y x R); thus condition (fl4"ll can be re-written as 
for all tp £ £>((-oo, 0) x R w ) smc/i Mai </J > 0, /or aZZ (t, x) £ (0, oo) x R N , we have 

M*<p(t,x) £ (Ci nc* 2 )°, 
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where C° denotes the normal cone of C. Let us recall that when the space dimension is finite (that is, if 
Ci, C2 are convex cones in R d for some d G N), then 

(cl(C 1 )ncl(C 2 )) =cl(C 1 + C 2 ), 

where cl(A) denotes the closure of the set A. 

If we forget about the closure and the fact that we are considering convex sets in an infinite dimensional 
space, then we are tempted to write 

M * ip(t, x) G (Ci n C 2 )° = Mi + A*2, 

with fa G C?, i = 1, 2. Moreover, very formally, we have 

C%"= "{d^m, m non-negative measure}. 

Thus, we may think of M. as some distribution of the form 

M = d^rn + fii, 

with m a non-negative measure on [0, 00) x M. xY xl, and /ii G C\. 

Of course, these computations are not rigorous, but we believe they may help the reader understanding 
the action of the distribution M. (at least in the divergence- free case), even though the precise structure 
of M. shall not be needed in the proof. Inequality lfT4|) is sufficient for all the applications in this paper. 

Let us stress that uniqueness for the limit problem holds, even though the cell problem does not have 
a unique solution in general; indeed, in the linear divergence free case, that is, if A(y,£) — a(y)£, with 
divj,a = 0, then a function u is a solution of the cell problem if 

div y (a(y)u(y)) = 0, (u) Y = 0. 

The constant function equal to zero is a solution of this equation, but in general there are other entropy 
solutions: think for instance of the case where N — 2, and 

0(2/1,2/2) = (-#2<Hyi, 2/2), di<Kyi,y 2 )), 

for some function cf> G Cp ei (Y). Then any function u of the form g{(f>) — {g{4>)) , with g a continuous 
function, is an entropy solution. Let us emphasize that nonlinearity assumptions on the flux are not 
enough to ensure uniqueness of solutions either, see for instance [18] . 

In Theorem [H the uniqueness of the solution of the limit system derives from a contraction principle 
associated with the macroscopic evolution equation, rather than the microscopic cell equation. The well- 
preparedness of the initial data, that is, the fact that uq(x, •) is an entropy solution of the cell problem, 
is fundamental. 

On the other hand, the lack of uniqueness of solutions of the cell problem entails that in general, 
there is no notion of homogenized problem. Indeed, if u is a solution of 

div y A(y,p + u(y)) = 0, (u) Y = Q, 

then in general, the quantity 

(A(.,p + u(-)) 

depends on u (except when N — 1, and in some special cases, when N — 2; see [13 [18]). Hence the 
macroscopic and microscopic scales cannot be decoupled: if l^< u (t,x,y) is a solution of the limit evolution 
problem, then u(t, x) = (u(t, x, •)) does not satisfy any remarkable equation. This is the main consequence 
of the absence of uniqueness for the cell problem. 

Let us mention an important particular case of the theorem [U which we call the "separate case". We 
now assume that the flux A can be written A(y,£) = ao(2/)s(£)> with div y ao = 0. This case has already 
been thorougly investigated by Weinan E in [9J in the case where g'(t;) ^ for all £, that is, when the 
function g is strictly monotonous. Here, we prove that his results hold with no restriction on g. 
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Let us introduce the so-called "constraint space" 

K := {/ £ L\Y)- div> /) = in &}, 

and the orthogonal projection Pq on K n L 2 (Y) for the scalar product in L 2 {Y). 
Then the following properties hold: for all f,g £ L 2 (Y), if / £ K , then 

Po(fg) = fPo(g). 

And if /, g £ Ko, then the product /g belongs to Ko. Notice also that all functions which do not depend 
on y belong to Ko. 

Proposition 1. Assume that uo £ L X (R N ,C per (Y)) PI L°°(M. N x Y), and uo is such that uo(x, ■) £ Ko 
for a.e. x E R N , 

Let 5,o = -Po(ao) £ L°°(Y). Let u — u(t, x; y) be the entropy solution of the scalar conservation law 

d t u(t, x; y) + div x (a (y)g(u(t, x; y))) = 0, (>0, i£ R N , y &Y, 
u(t = 0,x;y) = u (x,y). 

Then the function f(t,x,y,£) = l^ <u ( t x y } is the unique generalized kinetic solution of the limit 
problem lfT0|) with initial data l{< Uo (x,y)- In that case the distribution M. is given by 

M = ^ + 9 '(0(My)-My))-v x f, 

where m is the kinetic entropy defect measure associated with the function u, that is, f is a solution of 

d t f + g'(Oao(y)-Vxf = d i m. 
As a consequence, the solution u(t, x; y) of l|19p is an entropy solution of 

div y A(y, u) = 

for almost every (t,x) £ (0, oc) x R N . 
1.2 Convergence results 

Our first result is concerned with entropy solutions of (HJ . 

Theorem 2. Let A € W^ loc {R N+1 ). Assume that the initial data u a £ L] oc (R N ,C per (Y) satisfies lfT5]) . 
(fTBI , Let f = 1{< U be the unique generalized kinetic solution of the limit problem, with initial data lf< u ; 
the existence of f follows from theorem Ql Then as e vanishes, 

1-£<u E (t,x) ~ 1 l£<u(t,x,y)- (20) 

As a consequence, for all regularization kernels (p s of the form 
with ip £ V(R N ), J tp = 1, < <p < 1, we have, for all compact K C [0, oo) x R N , 



x £ R N , 



lim lim 

5->o e-»0 



u £ (t,x) -u* x if 5 (t,x,~) =0. (21) 



Remark 3. The assumption (fl5l) means thatuo is "well-prepared" in the sense thatuo(x, •) is an entropy 
solution of 

divy (A(y,u {x,y))) = 
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for a.e. x € K . // this hypothesis is not satisfied, then it is expected that the behavior of the sequence 
u e will depend on the nature of the flux. If the flux is linear, then oscillations will propagate, and the cell 
equation (9]) shall not be satisfied in general. If the flux satisfies some strong nonlinearity assumption, 
on the contrary, the conjecture is that the solution u £ re-prepares itself in order to match the microscopic 
profile dictated by the equation. Few results in this direction are known in the hyperbolic case; the reader 
may consult for instance JH, [791 [7^ \ 25]. In the author studies the same equation as (TJ) in which 
a viscosity term of order e is added, and proves such a result, but the method relies strongly on the 
parabolicity of the equation. 

Remark 4. The way in which theorem [H is stated might seem slightly peculiar; indeed, convergence 
results of the type 

u s — u (t, x,—^j — > in L] oc 
are expected to hold. But in order to establish such a result, it seems necessary to prove that 

lim I sup \u(t, x,y)-u * x y & (t, x, y) \dtdx = 0. 

But the evolution equation for u (or rather, for l£< u ) is given by definition^ since the distribution M. 
allows for very few computations, it seems difficult to derive such estimates. 

The next result generalizes theorem Q] to weaker solutions of equation (TJ, called kinetic solutions. In 
order to simplify the presentation, we explain how to generalize the result in the divergence-free case; we 
explain in the remark following the theorem how to derive an analogous result in the case where the flux 
A is arbitrary. 

Thus, for the reader's convenience, we first recall the definition of kinetic solutions in the divergence- 
free case (see [8] for the heterogeneous case, and the presentation in [23J for the homogeneous case) : 

Definition 2 (Kinetic solutions of (TJ)). Let u e G C([0, oo), L X {R N )). Assume that there exists a non- 
negative measure m e G C(M.£, w — M 1 ([0, oo) x R N ) such that for all T > 0, the function 

£ h- > / / rrf (t,x,^) dt dx 
Jo Jr n 

is bounded on M., and vanishes as |£| — > oo. 

Assume also that f e (t,x,t;) :— x(£, u e (t, x)) is a solution in the sense of distributions of the linear 
transport equation 



;)l , ^ Oi {-^)d Xi r = — t > 0, x e R N , (22) 
i=l ^ 

/ £ (i = 0)=x(e^o(^|)), (23) 

Then it is said that u e is a kinetic solution of equation (TJ). 

The existence of such solutions is only known when the flux satisfies additionnal regularity assump- 
tions. Assume that a, G Cp er (F x R) for 1 < i < N, and assume that there exists a constant C such 
that 

|o(i/,OI<C(l + |e|) VyerveeM. (24) 

Under such hypotheses, it is proved in [8] that for all uo G i 1 (M JV ,C pei (Y)), there exists a unique function 
u e e C([0, oo),L 1 (R N )) such that x{£,u £ ) is a solution of (TJ); u £ is called the kinetic solution of (TJ-©. 
And if u e is bounded in L°°((0, T) x R N ) for all T > 0, then u £ is the entropy solution of (TJ). Moreover, 
a contraction principle holds between kinetic solutions. 

Let us now state the convergence result for kinetic solutions : 
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Theorem 3. Let A G W*£? {oc (Y x R) such that div y A (y, £) = for ally,£. Assume that ^ G C\ er {Y xM) 
/or 1 < i < AT, and i/iai ([24)1 is satisfied. Assume that the initial data uq belongs to L 1 (Mr , Cp er (Y)) and 
satisfies 

N „ 

]T — (Oifo, Ox(£,«o))=0. 
i=l 

Lei u e € C([0, oo), L 1 (R Ar )) 6e i/ie kinetic solution of (JXJ) wii/i initial data Uo(x,x/e). Then there 
exists a function u G L°°([0, oo), L 1 (IR Ar x Y)) such that the convergence results lj20|) and l[2"Tj) /ioM, and 

Q^(a(y>Qx(€,u(t>x,y))) = in L> . 

Moreover, if we set 

d N d 

m := -x(e, u) + a *(y> O^rXU. «) G 

i—l 

then M. satisfies (fT8|) , 

Remark 5. Lei ws explain how this result can be generalized to the case where the flux A is arbitrary. 
First, the L 1 setting is not adapted to this case, because the L 1 norm is not conserved by the equation in 
general. Hence another notion of kinetic solutions is needed; the correct functional space should be of the 
type V + L 1 (R Ar ), where V is a fixed solution of the cell problem. 

Then, the crucial point in Theorem [3 is to find a sequence Uq such that Uq converges towards uq 
in L 1 (R Ar ,C P er(y)), and for all n E N, Uq satisfies ([HI , ((El). Finding such a sequence is easy in the 
divergence-free case, but seems more difficult in the general case, since solutions of the cell problem are 
not known. This seems to be the main obstacle to the generalization of Theorem^ to arbitrary fluxes. If 
this step is admitted, it is likely that the proof of Theorem [3 can be adapted to general settings. 

The plan of the paper is the following: in section [2] we prove, under the hypotheses of theorem [2j 
that the two-scale limit of the sequence l^ <u sf t>x \ is a generalized kinetic solution of the limit system. In 
section [31 we study the limit problem introduced in definition Q] and we prove the rigidity and uniqueness 
results in theorem [U hence theorem [Hand [2] will be proved by the end of section [31 In section 01 we study 
a relaxation model of BGK type, approaching the limit system in the divergence free case. In section 
we prove Proposition [TJ Eventually, in section [6l we have gathered further remarks on the notion of 
limit evolution problem. 



2 Asymptotic behavior of the sequence u £ 

In this section, we prove that the two-scale limit of the sequence f e — lg<u e (t,x)) sa Y / (t,%,y,£), is a 
generalized kinetic solution of the limit system; thus the existence result of Theorem Q] follows from this 
section. The organization is the following: we first derive some basic (microscopic) properties for the 
function f°. Then we explain how regularization by convolution can be used in two-scale problems. The 
two other subsections are devoted to the other properties of the limit system, namely condition lfl4ll and 
the strong continuity at time t = 0. 



2.1 Basic properties of /° 

We use the concept of two-scale convergence, formalized by G. Allaire after an idea of G. N'Guetseng 
(see [UEl]). The fundamental result in [l] can be generalized to the present setting as follows: 

Corollary 1. Let (g e ) 6 >o be a bounded sequence in L°°((0, oo) x M JV+1 ). Then there exists a function 
g0 g i°°((0,oo) x R N x Y x R), and a subsequence (e n ) such that e„ — > as n — > oo, such that 



g en (t,x,£)ip (t,x,-,n dtdx d£ ->/ / g°(t,x,y,£)'ip(t,x,y,£) dtdxdy d£ 

1 \ £ / ,/n ./B«vVvB 



R"+ 1 \ E / Jq JRNxYx 

for all functions if) G L 1 ((0, oo) x R N+1 ;C per (Y)). 

It is said that the sequence (g £ ")„ e N two-scale converges towards g . 
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Here, the sequence f £ is bounded by 1 in L°°\ hence we can extract a subsequence, still denoted by 
£, and find a function f° G L°°((0, oo) x R N x Y x R) such that (f £ ) two-scale converges to f°. It is 
easily checked that /° inherits the following properties from the sequence f £ 

0<f(t,x,y,O<l, (25) 
d^f° — —v(t,x,y,£), v non-negative measure. (26) 

Now, let us prove 0-©: let 

M := max(||ui||oo, IKIU) , 

where Ux,u% are the functions appearing in assumption l|16p . Since u$ (x/e) is a stationary solution of 
([1]), by a comparison principle for equation (TJ, we deduce that 

til (y — ^j < u £ (t,x) < U2 (y — ^j for almost every t > 0, x G 1 N . 

Thus ||ii e || L oo([ ,oo)xR JV < M, and for almost every t,x,£, for all e > 0, 

= i if£<-M, 

/ e (f,:r,£) = if£>M. 

Passing to the two-scale limit, we infer |(7j) and ((HJ). 

Now, we derive a microscopic equation for f°. First, multiplying |(4|) by S"(£), with S" G P(R), and 
integrating on (0,T) x B R x E, with T > 0, i? > 0, yields 

dt- 



a N+1 (-,t) f s S"(0dxd£dt = - [ [ [ m £ {t 7 x,^)S" (0 dx d£, dt, 

JRJB R V£ / Jo JRJB- 



where tir(x) is the outward-pointing normal to -B/? at a given point cc G 9 Br, and dctR(x) is the Lebesgue 
measure on SBr. 

Hence we obtain the following bound on m £ 

£to £ ((0,T) x B R xR)< C t ,r 

for all e > 0, R > 0, T > 0, and Supp m £ C (0, oo) x R x [-M, M]. 

Consequently, there exists a further subsequence, still denoted by e, and a non-negative measure 
m° = m°(t, x, y, £) such that em £ two-scale converges to m° (the concept of two-scale convergence can 
easily be generalized to measures; the arguments are the same as in [I], the only difference lies in the 
functional spaces). Moreover, Supp m°c(0,(X))xlxFx [—M,M], 

We now multiply (J4j) by test functions of the type eip (t,x,x/e,£), with tp G V peT ([0, oo) xR N x7xM), 
and we pass to the two-scale limit. We obtain, in the sense of distributions on (0, oo) x R N xYxR 

|rW^O/°)+|(aiV +1 ( y ,O/ )=^. (27) 

Thus is satisfied, which completes the derivation of the basic properties of f°. 
Now, we define the distribution 

+ f><*0f£ 

i— 1 

The distribution M. obviously satisfies ([6]). The next step is to prove that M. satisfies lfl4|): since regu- 
larizations by convolution are involved in condition (fl4|) . we now describe the links between convolution 
and two-scale convergence. 
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2.2 Regularization by convolution and two-scale convergence 

In this subsection, we wish to make a few remarks concerning the links between convolution and two-scale 
convergence. Indeed, it is a well-known fact that if a sequence (/„) weakly converges in L 2 (R N ) towards 
a function /, then for all convolution kernels tp — tp(x), the sequence (/„ * tp) two-scale converges in 
L 2 towards / * tp. It would be convenient to have an analogue property for two-scale limits. However, 
in general, if a sequence f £ — f e {x) is bounded in L 2 {R N ) and two-scale converges towards a function 
f = f(x,y) £ L 2 (R N x Y), then f e * tp does not two-scale converge towards / * x tp, Indeed, if %p = 
iP(x,y)eL 2 (R N ,C peI (Y)),then 



f e * <p(x)ip (x, —^j dx 

r(x'Mx-x')^ (x,-) 

v V e / 



dx dx' 



dx' f(x') 



(p(x — x')tp {x, — ^ 



dx 



In general, the quantity between brackets in the last integral cannot be written as a function of x' and 
x'/e, and it seems difficult to pass to the limit as e — > 0. 

In order to get round this difficulty, let us suggest the following construction, which is reminiscent 
of the doubling of variables in the papers of Kruzkhov, see |16^ I17j. With the same notations as above, 
consider the test function (tp * x tp) (x, f ), where <p(x) := ^p{—x) Va; G M. N . Then by definition of the 
two-scale convergence, 

f £ (x) [tp * x <p] (x, -) dx -> / f(x, y) [i> * x tp] (x, y) dx dy 



f e (x) [ip * x tp] (x, — ^ dx — J f £ (x')ip(x — x')ip ^x, — ~\ dx dx' , 



And 



f(x, y) [tp * x tp] (x, y) dxdy = / [/ * x tp] (x, y)ip{x, y) dx dy. 

R N xY JR N xY 

Consequently, as e — > 0, 

/ f e {x')tp(x - x')ip [x, — ) dx dx' -> / [/ * x tp] (x, y)ip(x, y) dx dy (28) 

JR 2N \ e J JR N xY 

for all tp e V(R N ), for all tp 6 L 2 (R N , C pei (Y)) . 

In fact, different assumptions on the function tp can be chosen; the key point is that tp should be 
an admissible test function in the sense of Allaire (see [l]). In particular, if there exist tpi £ V(R N ), 
tp2 £ L°°(Y) such that 

tp(x,y) = tpi(x)tp2(y), 
then tp is an admissible test function, and the limit (f28|) holds. 

2.3 Proof of the condition on M. 

The goal of this subsection is to prove that with 

N 

M = d t f + J2a i (y,0d i f, 

condition JW]) holds; hence, let tp £ P(R x E^), 6 £ V(R x R N ), such that 

tp > 0, 6 > 0, 

tp(t, x) = Vt > Vx E R N , 6(t, x) = 0Vt<0Vx£ R N ; 
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the function ip shall be used as a convolution kernel, and 9 as a test function, which explains the above 
hypotheses on the supports of ip and 9. 

Let ip € Q arbitrary (the definition of the set Q is given in definition [J) . We have to prove that the 
quantity 



oo poo 



A := I I I f(s,z,y,£) \d t ip(t-s,x-z) + y2a i (y,£)d i (p(t-s,x-z) 

10 JO JR 2Jv xYxK 



1=1 



x ^{Vi Q0(t> x ) d£, dy dx dz ds dt 



in non-positive. 

Before going into the technicalities, let us explain formally why the property is true; let us forget 
about the convolution and the regularity issues, and take the test function 



in equation Q. 

Let R > max(M, C); recall that M and C are such that Supp f° C [0, oo) x R N xY x [-M, M), and 



ip(y,£) — a- if £ < — C. Integrating on [0, oo) 



piV 



[-R, R], we obtain 



o Js. N J-R 

oo r pR 



1 

£ Jo 



f E (t,x,0 [d t 6{t,x)+a i (^,£)d Xi 0{t,x)\ ^(^C) dxd£ 

f s (t,x,o^(^,t)e(t,x) dxd^dt 



o Jr n £ 

OO p rR 

/ / m 

JR N J-R 



-a^v+i ^ — , — Rj 0(t, x) dt dx 
(s, z, £)d£ip ^— , £j dz dt; ds 



£<u (a;,f ) 



9(t = 0,x)ip dxd£. 



Notice that 



^a N+1 -Rj = -div x A -R 



and thus 



JR N J-R 
oo p pR. 



f E (t,x,0 [dt6{t,x)+ ai (~,^)d Xi e(t,x)\ dxd£dt 



(p -R) 3,9(1, x) dt dx - J J l e<Uo ( x , f ) 0(t = 0, x)ip ( J, e) ^ ^ 

/ / A i (-,-R)d i 9{t,x)dtdx- ( [ l* <uo ( a *)0(t = 0,x)ip(~,t) dxdi. 
Jo Jw N ye 7 Jr n J-r v ' e; Ve y 



x) <i£ ds 



-a. 



-n / / A 

'0 Jl» Ve 



Passing to the limit as e — > 0, we retrieve 

,-R 



f°(t, x, y, [$0(f, x) + a, (y, ^(t, a:)] V (», 6 dx dy d£ dt 

JR N J-R. 

-n / / A l (y,~R)d l 9(t,x) dtdx - / / l 5<Uo(:Eia) 0(i = 0, x)V> (y, d£ 
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This means exactly that 



d_ 

dt 



Oi/V < o, 



or in other words, that J Yx r < in the sense of distributions on [0, oc) 



Now, we go back to the regularizations by convolution. According to the preceding subsection, 



A = lim 



OO fOO 



JO 



f e (s,z,0 \ d t <p(t -s,x-z) + S2ai di<p(t - s,x- z) } x 

x ^ (~>^) ^(*> x ) d£ dx dz ds dt. 

Hence, in |(4]), we consider the test function 



(j)(s,z,€) 



JR N 



ip(t — s,x — z) 6(t, x) dt dx 



where K is a cut-off function such that < K < 1, K £ = 1 if |£| < R, and 

V\5 := 1p *y if{ *£ <^2) 

with e PQR^), ^2 £ £>(»), < y>j < 1, J" <ft = 1 for i = 1, 2, and 



rffo) - ^ (f) , m = ^ (| 



According to (j4]) , we have 



JK»+! 

I /"°° 

^ Jo Jr n + 



N 

d s (t>(s, z, £) + ^2 a * (^) <9 Z ,<K S ^£) 



ciz <i£ (is 



JR N + 



i / e (s,z,O a JV+i d £<t>{ s ,Z,€) dzd^ds 

m e (s, z,t;)d£(j)(s, z,{;) dz d£ ds + / X k, ( z, — ) J 0(s = 0, z, £) dz d£ 
i itiv+i V v eJJ 



(29) 



And 



d s <f>(s,z,g) = 
V z 0(s,z,£) = 



Jm. n 



+- 



JK" 



dtip(t — s,x — z) 6(t, x) dt dx 
\7 x ip(t — s.x — z) 9{t, x) dt dx 
ip(t — s, x — z) 9(t, x) dt dx 

N 

tp(t — s, x — z) 6(t, x) dt dx 
ip(t — s,x — z) 0(t, x) dt dx 

4>s K{£) = 0. 



JR N 

dtf(t, x — z) 0(t, x) dt dx 



^ (f,e)if(o, 

(v v ^)(pf)if(0, 



/o Jk n 

Thanks to the assumption on the support of tp, and the fact that 

dfips = (d^ip) * y <Px *£ (p s 2 > 0, 
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we have 



JR N 



ip{t — s, x — z) 9(t, x) dt dx 



>0. 



Moreover, thanks to (7]), ([8]), and the assumptions on ip and K, we have d^K = on Supp m e , and 



o Js. N 



ip(t — s,x — z) 0(t, x) dt dx 



la Jr n 

Hence, we obtain, for all e,S > 0, 



ip(t — s, x — z) 6(t, x) dt dx 



WO- 



-J f £ (s, z i£) ^d t tp{t - s,x - z) + Y2 ai (f'O ditp ^ ~ s ' x ~ z )| x 

xif>S (^ — , tij 0(t, x) d£ dx dz ds dt 
+ f £ {s,z,£) a -Vy^S (~>f) f{t- s,x- z)9(t,x) K(^)dtdxdsdzd£_ 

H / </?(f — s, x — 2) 0(t, x) 9^i"T(^)aAf + i ^ — , £j dt dx ds dz dt; 



0. 



Following the formal calculations above, we have to investigate the sign of the term 

f £ (s, z, £) a ' ^y.^s ^{t — s,x — z) 6{t, x) K(£) dt dx ds dz d£. 



Since div y ^(aip) = —d^/i, we have 



dWy^{aip s ) = -~7jr + r s 



where fis — fi * y ip{ ip\ ■ Then 

.ft r ^o a -§(-,()\r I 

Jo ot, \e / [Jo Jr 

^ = ii e (t,i)) W (^) 



<^(t — s, a; — z) 0(4, x) di dx 



dz d£ 



JR N 



tp(t — s,x — z) 6{t, x) dt dx 



ds dz d£ < 0. 



'0 JR N + 1 

Hence, we have to prove that as S — > 0, 

r s -> in Lj^F x R). 

The proof is quite classical. We have 

r 5 (y, = a(y, £)V> * (VyM<P%) ~ Wv> £M * (V tf> «¥>M) 

iV 

= X] / M^C) ~ a i(yi'Ci)] 1>(vi,€i)dyi<pi(v - yi)A(^ - 6) ^1 

2—1 

+ / [«iv+i(y,0 -ajv+i(yi, Ci)]'0(yi, - yi)5«^£ - £1) */i *a 
Thus, we compute, for G JR 2A,+2 , 1 < i < JV + 1, 

- = (y-Vi)- W v ai{ry + (1 -t)j/i,t£ + (1 -t)£i) dr 



+ (C-a)' / ^(Ty+Cl-rJjtt.rC+Cl-T^iJdT. 
Jo 
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Set, for 1 < k, i < N, y € R N , £ G 



Notice that 



</>fc,i = — 5fe,i, / Ci — Sff+i^. 

f+i Jb jv + 1 

Then 

JV+l iv . „ 

r«(i/,0 = EE/ + + (1-^)6 Wl/i,Ci)4i(»-I/i,C-ei)di/ideidT 

+ Ey ^(ri/ + (l-r)yi,re + (l-r)e 1 )^( yi> $ 1 )^(»-l/i,e-6)dl/ideidT. 
Hence as (5 — > 0, converges to 

-div l/ , e (o(y,0)^(»,0=0 

in L|' oc (R Ar+1 ) for any p < oo and for all (i, x) £ [0, oo) x R N . We now pass to the limit as 6 — ► 0, with 
e fixed, and we obtain 

/ e (s, z, £) |e?{<£(£ — s, a; — z) + ai ^— , £^ 9j<^(t — s, a; — z)| ^ ^— , £j 9(t, x) d£ dx dz ds dt 

J 9(t,x)d i K(^)A^-,^j -V x (p(t - s,x - z) dtdx ds dzd£ 



> 0. 

Passing to the limit as e vanishes, we are led to 

/°(s, z, y, f) {d t (p(t - s,x - z) + di (y, £) di<p(t - s, x - z)} tp(y, £)9(t, x) d£ dx dz ds dy dt 
-a_ I 9{t, x) d^K(^)A (y, £) • V x <p(t — s,x — z) dt dx ds dy dz d$, 



> 0. 

Since 

J Hi, ,)W ---=)<* * * * ~ (/ *. *) * *) (/ v.,(», ,)**)- o, 

we deduce that 

J f{s,z, y,£) - s,x - z) + E a » - s,a; - ^)| ^ (y,£) 9(t,x)d£ dx dz ds dy dt < 0, 

which means that /° satisfies condition l[T4"ll . There only remains to check the strong continuity of / at 
time t = 0. 

2.4 Strong continuity at time t — 

The continuity property for/ is inherited from uniform continuity properties at time t = for the 
sequence f £ . This is strongly linked to the well-preparedness of the initial data (condition (J9j) ) , that is, 
the fact that for all x £ M. N , uq(x, ■) is an entropy solution of the cell problem 

div y A(y,u (x,y)) = 0. 
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First, let us consider a regularization of the initial data 



9n = fa *x Pn *y <P\ *£ <£>2- 

with p n a convolution kernel (n G N), 6 > 0, and tpf defined as in the previous subsection. Then we can 
write 



i=l 

1 /a; 
-a - 

e Ve 



(30) 



Notice that 
and 
where 



\Vx9i\\L°°(g. N xYxM) < I|VxPti||l1 

a (2/, ^ydi {x, V, = <% m n + r*, 



= m Q P„ *J, *£ (^2, 



r n {x,y,i) =a{y,i) V y g n {x,y,£) - [af„ * x p n ] * y>£ V^c^y)^^)- 

Then for all neN, for all a; G R N , r s n vanishes as 5 — > in Lj oc (Y x R) and almost everywhere. The 
proof of this fact is exactly the same as in the preceding subsection, and thus, we leave the details to the 
reader. As a consequence, 

and there exists a constant C n , independent of e, such that for all n G N, for all e > 0, and for almost 
every x, £ 

limsup|i£ { (a;,£)| < C„. 

Moreover, Supp R e n S C R N x [-i? - 1, R + 1] if (5 < 1. 

Now, we multiply I® by 1 - 2g n (x,x/e,£), and ([3U]l bv 1 - 2f £ {t,x,£). Setting 

/i*, 5 (t,z,0 := / e (*,*,0 [l - 2g 5 n (x,^,^)] +g s n (x,^,t) [l-2f(t,x,0} 

= f e (t,x,0-gi(x,^,t)\\gi(x,^t)-\g s n (x,^z) 

we obtain 
d 



dt 



-hi 



N 1 x 

i—1 



dm € 
Notice that 



1 - 2 5 £ (ar, p e)] + ~%m* (x, V s {t, x, 0] + K,s( x > I 1 " 2 / £ ( f > 0] • (31) 



a ? [l-2/ £ (i,x,C)] =2<5(e = u £ (t,a;)), 
4(l-2<£ (*,^))=2,„, e , aM , 

where v n ,e,s is a non-negative function in C°°(R Ar+1 ), with support in M. N x [-M — 1,M + 1] if S < 1. 
Notice also that f e (t, x, £) — g n (x, x/e, £) — if |£| is large enough (|£| > A/ + 1). Take a cut-off function 
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C = ((x) such that C(x) = e 1*1 when |x| > 1, and i < £(x) < 1 for \x\ < 1. Then there exists a constant 
C such that 

|VxC(*)l <CC(x) Viel" 
Hence, mutliplying ((3TJ) by £(x) and integrating on we obtain a bound of the type 

d_ 

lit 



K,s(t, %,QC(z) dx d£ < C h* s (t,x,£)t(x) dx d£ 



+ R e n , 5 (x,0 \l~2f(t,x,0\C(x)dxdt. 

Using Gronwall's lemma and passing to the limit as S — > with e and n £ N fixed, we retrieve, for all 
t > 0, 



RN + l 



f e (t,x,£) -g n (x, ({x)dxd£ 



< e 



ci 



,ci 



So (x, -,£) -g n (x, 



(x, - 5n (a;, j,£J 



£(x) d£ 



£(x) dx e?£ 



+ C n (e ct -1), 

where the constant C n does not depend on e, and g n — /o * x Pn- And for all n £ N, e > 0, we have 

So (x, ~, - g n (x, C(x)dxd£ 

So (x, - So (x 1 , p n (x - x')C(x) dx dx' d£ 

"0 (xi--;£) ~ u Pn{x ~ x'%(x) dx dx' 



< 



< 



< 



sup \uo (x, y, £) - u (x', y,0\Pn(x - x')((x) dx dx'. 

R N y£Y 



RiV + 1 



The right-hand side of the above inequality vanishes as n — > oo because uq £ Lj 0C (M iv , C pex (Y)). Similarly. 

9n (x, - g n (x, j,^ ((x)dxd£ 
g n (x, - So (x, ((x)dxd$, 

5n (z, j,£J - /o j,^) C(x)dxd£, 
sup |u (x, y, - wo (a;', y, £)| P«(x - x')C(z) dx rf a;'. 



< 

< 3 
3 



Hence, we deduce that there exists a function w : [0, oo) — ► [0, oo), independent of e and satisfying 
lim t ^o = 0, such that 



RN + 1 



RN+1 



S e (t, x, - /o (x, p C) I C(x) rfa; d£ < w(t) 



for all t > 0. 

Then, we prove that the same property holds for the function / . Indeed, we write 



S e (t,x,0-i 



£<u (x,f ) 



f - 2fl 



?<«o(s,f ) + £<wo(a:,f )' 
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let 6 £ L°°([0, oo)) with compact support and such that > 0. Then for all e > 0, 

^ [r 2ri i<uo{xf) + i s<uo(x , f )] ewe® dx dt dt < ^ Lo(t)9{t) dt. 

Since uq G L 1 1 oc (R Ar , C pei (Y)), it is an admissible test function in the sense of G. Allaire (see [1]); we 
deduce that 1^< U0 is a l s ° an admissible test function. This is not entirely obvious because it is a 
discontinuous function of uq. However, this difficulty can be overcome thanks to an argument similar to 
the one developed below in subsection 13.31 and which we do not reproduce here. Thus, we can pass to 
the two-scale limit in the above inequality. We obtain 

(/"(*, x, y, - \f(t, x, y, 0| 2 + \f(t, x, y, £) - l (<Uo{x , y) \ 2 e(t)((x) dt dx dy d£ < 

/■oo 

< / 9(t)u(t) dt 
Jo 

Notice that |/°| — \f°\ 2 > almost everywhere. As a consequence, taking 9(t) = lo<t< T , with r > 
arbitrary, we deduce that 

-/ \f(t)- X (^uo(x, y ))\ 2 C(x)dtdxdy<- f oj(t)dt, 
T Jo T Jo 

and the left-hand side vanishes as r — > 0. Thus the continuity property is satisfied at time t = 0. 

Hence, we have proved that any two-scale limit of the sequence f £ is a solution of the limit system. 
Thus the existence result in Theorem [TJ is proved, as well as the convergence result of Theorem [2j We 
now tackle the proof of the uniqueness and rigidity results of Theorem [TJ The strong convergence result 
of Theorem [TJ will follow from the rigidity. 



3 Uniqueness of solutions of the limit evolution problem 

In this section, we prove the second and the third point in Theorem [TJ that is, if / is any solution of the 
limit evolution problem, then there exists a function u G L°°([0, oo) x R w x Y) such that f(t,x,y,£) = 
l£<u(t,x,y) almost everywhere, and if f± = lf< Ul , /2 = l£<u 2 are two generalized kinetic solutions, then 
the contraction principle (fT7|) holds. 



3.1 The rigidity result 

Let / be a generalized kinetic solution of the limit problem, with initial data 1j<m . The rigidity result 
relies on the comparison between / and f 2 . Precisely, we prove that / = |/| 2 almost everywhere, and 
since <%/ = — v < 0, this identity entails that there exists a function u such that / = lg<«. Thus, we 
now turn to the derivation of the equality |/| = |/| 2 . 

Let 8 > arbitrary, and let Q x G X>(R), 6 2 G V(R N ) such that 

0i > 0, 9 2 > 0, 

l 0i = l 0a = 1, 
Supp(9i c [-1,0] and 6>i(0) = 0. 

We set, for (t, as) £ R N+1 




Set f := f * t ,x d S , M s := M * t ,x s . Then f s is a solution of 

l — l 
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Moreover, f s satisfies the following properties 



< f < 1, (32) 
liv y4 (a(y,£)f s )=d i :m* t , x d 5 , (33) 
d s f = -v * t , x e 5 , (34) 



whereas A4 satisfies 



M s G C((0, T) x E w , L 2 (F x R)) n L°°([0, oo) x x Y x R), (35) 

M s (-,0=0 if ICI >M,f(-,O=0 if£>Af, = l iU<~M, (36) 

M s ip<0 VipeG. (37) 

K 

In particular, notice that (1 - 2f s (t, x)) G Q for all t, x, and f 5 (t, x, y, f) - / 5 (t, x, y, £) 2 = if |£| > M. 

Let C G C°°(IR Ar ) be a cut-off function as in the previous subsection. We multiply the equation on / 
by (1 — 2f s )£(x), and we integrate over R N xFxR, We obtain 



| / (f s -\f\ 2 )c- f at(yma*){?-\f\ a )= I 

al JR N xYxM JR N xYxR JR 1 



M (l-2/ d )C <0. 

xYxR 



We then deduce successively, using Gronwall's lemma, 



d 
dt 



J 

Jr.' 



(f s -\f\ 2 )C<c (f-\f 5 \ 2 )C 

'xYxR JR N xYxR 



(f(t) - \f(t)\ 2 ) (<e ct {f(t = 0) - \f(t = 0)| 2 ) C V<> 0, 

<YxR JR N xYxR 

T r p ct - i r 

/ (f 5 \f 5 \ 2 ) C < / (f(t = 0) - \f(t = o)| 2 ) C, (38) 

VxYxR ° JR N x Y xR 

with a constant C depending only on ||a||ioo(y x r_ 

Let us now check that f s (t = 0) strongly converges towards 1^< U0 = /o at time t = 0. In fact, 
the main difference between the proof of Theorem Q] and the one for generalized kinetic solutions of 
scalar conservation laws (see chapter 4 in [23] ) lies in this particular point. Indeed, in the case of scalar 
conservation laws, the continuity property can be inferred from the equation itself; in the present case, 
the lack of structure of the right-hand side M. prevents us from deriving such a result, and hence the 
continuity of solutions at time t — is a necessary assumption in definition [TJ 

Using hypothesis lfl3|). we write, for almost every x,y,£, 

f s (t = 0,x,y,0 = f f( Sl z iy ^)6 s (~ Sl x~z)dsdz 



<r+i 



f(t = 0,x,y,S)-f * x e? 2 (x,y,0 = / (f(s, z,y,£) - f (z,y, 0) d (-s, x - z) ds dz. 

JR N + 1 

As a consequence, for all S > 

\f 5 (t = 0)- f * x 6 s 2 \ 2 ((x) dx dydH 



j 

.75 



< / / \f(s,z,y,£) - fo(z,V,Q\ C{x)0 {~s,x~ z)dxdyd£,dsdz 

JR N xYxR JR« + i 

< ll/(*)-/o|li»(R Wx yxR,c(x)«ix < iv«if)^i (x) dsdxdydt + 2R\Y\ ||C-C*^II^ 
C f s 

~ ~S Jo \\ f ^~ f0 ^ 2 ^ N xYxR,C(x)d xd ydO ds + 2R \ Y \ IIC-C*02llz*(B*)' 
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The right-hand side of the last inequality vanishes as 5 — > 0, and thus f s (t = 0) converges towards /o as 
5 — » in L 2 (R JV x7xl, £(x) da: dyd£), and hence also in L 1 (M Ar x7xl, ((x) dx dyd£). Consequently, 

(/ 5 (* = 0)-|/ 5 (t = 0)| 2 K^0 as^O. 

Above, we have used the fact that fo — 1^<« , and thus fa = f§. 

Now, we pass to the limit as S — > in (|38|) : we obtain, for all T > 0, 

(f-\f\*)<P<o. 

<YxR 

Since the integrand in the left-hand side is non-negative, we deduce that |/| = |/| 2 almost everywhere. 
The rigidity property follows. 

3.2 Contraction principle 

Let /i, /j be two generalized kinetic solutions of the limit problem; we denote by Mi, M2, and .Mi, .M2, 
the constants and distributions associated to /1, /2, respectively. Without loss of generality, we assume 
that Mi < M 2 . According to the rigidity result, there exist functions ui,u 2 G L°°([0, 00) x R N x Y) n 
i oo ([0,oo),L 1 (R Ar x Y)) such that /; = l c<tlt . 

As in the previous subsection, we regularize /i, A^i by convolution in the variables t, x, and we denote 
by ff, Mf the functions thus obtained. The strategy of the proof is the same as in [23], Theorem 4.3.1. 
The idea is to derive an inequality of the type 

- J \f 1 (t,x,y,Z)-f 2 (t,x,y,0\({x)dxdydfi<C J \fi(t,x,y,£) - f 2 (t,x,y,£)\C(x) dxdy d£, (39) 

where £ is a cut-off function as in the previous section. 

Since |/i(£) — f2(t)\ = \fi(t) — y 2 ( ^ ) | 2 = fi + fi — 2/1/2, let us first write the equation satisfied by 
9 5 ■■= ft + fi - Vifl We compute 

\d t f( +^ ai (y,£)Aff=.M?j x 1-2/1, 



i=i 



Adding the two equations thus obtained leads to 

N „ 

dt9 5 + E °*(tf> O^V = [1 - 2 /l] + M [1 - 2/f] • 

i— 1 

Notice that thanks to (0, iJH) and the microscopic constraints ([9]), (fTTjl . 1 — 2ff(t,x) E 5 for all (i, x). 
Hence 

A4|(t,aj) [1-2/f (*,»)] <0 V(t,x) E [0,oo) x R N , 

R 

and the same inequality holds if the roles of f\ and /2 are exchanged. 

Now, take a cut-off function C £ C°°(R ) satisfying the same assumptions as in the previous subsec- 
tion; multiply the equation on g s by ((x), and integrate over R N x Y x E; this yields 



d 
dt 



and thus 



/ 



9 S (t, x, y, £)C(x) dxdyd£<C / g 5 (t, x, y, £)C(x) dx dy d£ Vt > 0, 

'xFxK JR N xYxR 



g s (t, x, y, 0C(z) dx dy d£ < e ct / g s (t = 0, x, y, £)C(x) dx dy d£. 

xYxR JR N xYxM 
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According to the strong convergence results of ff(t = 0) derived in the previous section, we can pass to 
the limit as S — > 0. We infer that for almost every t > 0, 

/ \fi(t,x,V,0 - f2(t,x,y,£)\((x) dxdydS, 

< e ct f \fi(t = 0,x,y,Z)-f 2 (fi = 0,x,v,Z)\C(x)dxdyde. (40) 

JM N xYxM 

This completes the derivation of the contraction principle for the limit system. Uniqueness of solutions 
of the limit system follows. In particular, we deduce that the whole sequence f £ of solutions of Q two- 
scale converges towards / . 



3.3 Strong convergence result 

Here, we explain why the strong convergence result stated in Theorem [2] holds, that is, we prove l|2ip . 
This fact is rather classical, and is a direct consequence of the fact that 



2 sc. 



*-£<u(t 1 x 7 y) • 



Let us express this result in terms of Young measures: the above two-scale convergence is strictly equiv- 
alent to the fact that the two-scale Young measure Vt,x,y associated with the sequence u £ is the Dirac 
mass <5(£ = u(t, x, y)) (see [23], Chapter 2). And it is well-known (see [9]) that if u is a smooth function, 
then 



dv t , x ,y{0 = 6(£ = u(t,x,y)) 



u — u 



(t,x,^j^0 in L\ oc . 



For the reader's convenience, we now prove the result without using two-scale Young measures. We 
define us — u * x (p§ , with ip$ a standard mollifier. Take K G T>(R) such that < K < 1, and K(£) = 1 if 
|£| < M. Take also a sequence Q n G C°°(R) such that < 6 n < 1, and 

e n {i) = !)!(,<--, e n (o = oifc>i. 

n n 



Then we write 



^ min(u e (t ,x) ,u$ (t.x, &\\<£<xnax(u e (t,x),U8 (t,x, j)) ' 



The function 1 
replace it by 



£<u e (t,x,f) 



is not smooth enough to be used as an oscillating test function. Thus we 



?n (S, - US (t,X, ^jj 



and we evaluate the difference : for all compact set C C [0, oo) x 

h<u s (t,x,f) " °- " u * (*> *> f)) | K & dt dx dC <l\C\. 
According to the two-scale convergence result, for all n G N, 



£-u s [t,x, -J J li <U '(t,x) K (0 dt dx d£ 



n (£ - u s (t, x, y)) l£ <u ( tiXty )K((;) dt dx dy d£. 

Since the sequence 9 n (£ — us) uniformly converges towards l^< Uj as n — > oo, we can pass to the limit as 
n — > oo, and we deduce 

1 £<u ~(tx ^)h<w(t,x)K(^) dtdxd£ -> / / l i : <Us{t , x ^ ) l i : <u{t . x . y) K{C)dtdxdyd£ > . 
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Simlarly, as e — > 0, for all 6 > 0, 

/ / 1 £<uJtx^) K (O dtdxd C^ / / lz<u s (t.x.y)K(£)dtdxdyd£„ 
Jc Jr. v ' ' " ' Jc Jrxy 

/ / lt<v*{t,x)K(£)dtdxd£-> l s<u ( tiX>y) K((,)dtdxdyd£. 
JcJrxy Jcjr 

Thus 



K(£)dtdxd€^y / / | le<«(t,x, y ) - !£<u«(mmoI rf * dxdyd£ 

Jc JRxY 



On the other hand, 



L?<««(t,x) - 1 4 <„,( t ^,|) K(£)dtdxd£= u E (t,x)-u s \t,x,-J 



Li(C) 



c Jixy 



*-£<u(t,x,y) 



lz<u s (t,x,y)\ K(0 dtdx dyd£=\\u- us\\ L i(CxY)- 



Hence we have proved that for all 5 > 0, for all compact set C C [0, oo) 



pJV 



lim 



u e (t, x) — Us (t, X, — J 



\ u - usWl^cxY)- 



Statement ll2"Dl then follows from standard convolution results. 



3.4 Application: proof of the convergence result for kinetic solutions 

In this subsection, we prove Theorem [3J this result is in fact an easy consequence of the convergence 
result stated in Theorem [2] for entropy solutions, and of the contraction principle for the limit system. 
Assume that ajv+i = 0, and let u e be a kinetic solution of equation ([1]), with an initial data uq{x,x/s) 
such that u a G L l (R N , C pei {Y)) and 







TT- & x (€> y ^ = 



in the sense of distributions. 

For neN, let uft := sgn(u ) 'm£(\u \,n). Then for all n G N, wg belongs to L°°(R N x Y) and 



(41) 



uq as n — > 00 in L 



CperPO)- 



Moreover, x(£j u o) = x(£, Uo)l(;<nj and thus for all n G N, Uq satisfies (|4ip . 

For all n,e > 0, let < G C([0, 00), L 1 (R iV )) n L°°([0,oo) x 1*) be the unique entropy solution of 
equation (JXJ with initial data Uq (x, x/e). Then by the contraction principle for kinetic solutions of scalar 
conservation laws, we have 

Vn G N, |K - <||l~([o,oo),li(k«)) < u (x, - Uq (x, ^ ^ < ||u - Uo IU^rw,^^))- 

On the other hand, for all neN, let If < u „ be the unique solution of the limit system with initial 
data l$<u%- By the contraction principle for solutions of the limit system (see inequality (|4"0)l ). we have, 
for all integers n, m G N, for all i > 0, 



|u„(t,a;,y) - u m (t,a;,y)| f(x) dx dy < 



'xY 



W™{t,x,y) -Uo(t,x,y)\((x) dx dy 



xY 



cC*ll„, m ,,"11 

< e ||u -u || L i( K w iCper (r)), 



(42) 



where C G C°°(R JV ) is a cut-off function satisfying the same hypotheses as in the previous subsections. 
Consequently, the sequence (u n ) n ^ is a Cauchy sequence in L^ c ([0, 00), L l (R N x Y)); thus there exists a 
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function u G Lj^ c ([0, oo) , L 1 (R N xY)) such that u n converges towards u as n — > oo in ij^ c ([0, oo), L 1 (R W x 
Y)). Moreover, the limit u is independent of the chosen sequence Uq thanks to f42|) : indeed, if Vq , Wq 
are two approximating sequences giving rise to functions v and w respectively, we construct the sequence 



v 







»o if n is even, 
Uq if n is odd. 



Then the sequence Uq converges towards uq, and thus the corresponding sequence u n converges towards 
u, while U2n converges towards v and U2 n +i towards w. By uniqueness of the limit, u = v = w. 

On the other hand, since the sequence f e — x(£, u e ) is bounded in L°°, there exists a sequence (ek)k&i 
of positive numbers, £k — ► 0, and a function / 6 L°°([0, oo x x Y x R), such that 

x(£, (*, x)) /(*, x, y, f). 

Now, for all fe, n £ N, 

||x(£,u e *) ~ x(C< !c )IIl~([o,oo),l 1 (r«+ 1 )) < IK - u ollL 1 (R JV ,c per (y)) I 
and for all n G N, since x(£i w ) — lf<« — l?<o, we have, as k — > oo, 

Let G X> per ([Q, oo) x R w x Y x R). By definition of two-scale convergence, 
[ X (£, w £fc (t, x)) - < fc (t, a;))] y> (i, a;, — ,£] dt dx d£ -> 

/(*> ac, y> - x(£= u n(*> x > v))\ V (*> 2/> ^ d 2/ d £- 

o ii"xyxi 
And for all fe G N, the following inequality holds 

X \ 

< 



f [x((,« £i (M))-x(e,< k (t,i))]Ht,i,- ,e) d^xdc 

< H'^llL 1 ([0,oo),L°=(R N xyxR))ll u O ~ u I |i 1 (K 7 ',Cp or (y)) ■ 

Passing to the limit as k — ► oo, we deduce that for all n G N, <p G P pe r([0, oo) x R N xYx R), 



[f(t, x, y, £) - x(£, «n(*> y))] 93 (t, y, £) di da; d£ 

'xYxR 



< 



< |K — WolU 1 (K JV ,C psc (> r ))ll¥'l|L 1 ([0,oo),i«>(R N 'xrxR)- 

Thus, we pass to the limit as n — > 00 and we infer that / = x(£, M (^ a;, y)) almost everywhere. Hence the 
limit is unique, and the whole sequence x(£, u e ) converges (in the sense of two-scale convergence). 

Eventually, let us pass to the limit as n — > 00 in the limit evolution problem for x(^ u n)- We set 
/ = x(£j m )> an d define the distribution 

M :=d t f + a(y,0-V x f. 

Then M n — M. in the sense of distributions, and it is easily checked that inequality lfl8|) is preserved 
when passing to the (weak) limit. Thus A4 satisfies lfT8|). 

In the divergence-free case, the main difference between the L°° and the L 1 setting, that is, Theorem[2] 
and Theorem [3l lies in the fact that uniqueness for the limit system in the L 1 setting seems difficult to 
derive; indeed, the proof of uniqueness in the L°° case uses several times the fact that the distribution 
M has compact support. In a L 1 setting, this assumption would have to be replaced by a hypothesis 
expressing that A4 vanishes as |£| — » 00, in some sense. But it is unclear how to retrieve such a 
property from the hydrodynamic limit (see section 2]), for instance. The above argument only proves 
that uniqueness holds among L 1 solutions which are obtained as the limit of a sequence of L°° solutions. 
Thus we have left open the correct notion of limit system in a weak L 1 setting, and the derivation of 
uniqueness therein. 

Nonetheless, we wish to stress that the contraction principle in the L°° setting is sufficient to ensure 
that the whole sequence x(£; u£ ) converges, even if uniqueness for the limit system fails. 
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4 A relaxation model for the limit evolution problem 



In this section, we exhibit another way of finding solutions of the limit system in the divergence-free case. 
Indeed, the existence result in theorem [TJ was proved by passing to the two-scale limit in |(4|), and it may 
be interesting to have another way of constructing solutions, which does not involve a homogenization 
process. 

Hence, we introduce a relaxation model of BGK type, in which we pass to the limit as the relaxation 
parameter goes to infinity. The drawback of this method is that the existence of solutions of the limit 
system is not a consequence of the construction. Indeed, we shall prove that if a solution of the limit 
system exists, then the family of solutions of the relaxation model strongly converge towards it in the 
hydrodynamic limit. Hence the proof is not self-contained, because the existence of a solution of the 
limit system is required in order to pass to the limit. Nevertheless, the final result may be useful in other 
applications. 

In the whole section, when we refer to the limit system, we have in mind the modified equations 
introduced in Remark [TJ In the divergence-free case, it is also slightly more convenient to work with the 
function x(£, u), rather than lg<«. Hence a solution of the limit problem is a function g satisfying 

N ~ 

£~-My,a9)=0, (43) 

d i9 = 6(0-v(t,x,y,0, ^>0, (45) 

and M is such that for all tp G P([0, oo) x R N ) such that ip > 0, the function M * t ,x <P belongs to 
C([0, oo) x M w , L 2 (Y x K)), and 

(-M *M ^) (M, •) V> < 0, 

' YxR (46) 

VV> G L£ C (Y x M), div y (aip) = 0, and d t ip > 0. 
4.1 A relaxation model 

The goal of this subsection is to introduce a system approaching l(43|) - (|46l) . With this aim in view, we 
define a relaxation model of BGK type, which takes into account the constraints the limit system, that 
is, equations (|43 ]l -(|46" ll . Let 

M := ||wo||l~(yxr), 
E := {/ G L 2 (Y x M), Supp/ C Y x [— M, M]}, 
K := {ip G E, div y (a(y, t)<p(y, 0) = in V .} , 
JC := K n {tp e E, Bu G M^ et (Y xR), u> 0, d^ip = 5(0 - v) . 

Then E endowed with the usual scalar product is a Hilbert space, and K is a closed convex set in E. 
Thus the projection V on fC is well-defined. 

The main result of this subsection is the following : 

Proposition 2. Let X,T > be arbitrary. Set 

X T :=C([0,T],L 2 (R^ xYxR s )). 

Then there exists a unique solution f\ G Xt of the equation 

dth + a(y, • V x f x + A/a = XV(f x ), (47) 
h(t = 0,x,y,£) =x{€,uo(x,y)) 

The function fx has the following properties : 
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1. For almost every t,x,y,^, 



fx{t,x,y,Z) = OifS>M, 
sgn(Of\(t,x,y,0 = \ f\(t,x,y,g)\ < 1. 

2. L 2 estimate: for all A > 0, 

\\f\\\x T < H«ollLi(R«xy)- (48) 

3. Strong continuity at time t = 0: there exists a function u> : [0, oc) — > [0, oo), such that lim + u> = 0, 
and such that for all A > 0, t > 0, 

1 1 AW - /ollL'(R.vxyxR) < w(t). (49) 
4- Fundamental inequality for A4\ := X(P(f\) — f\): for all g £ IC, for almost every (i, x), 

M x (P(h)-g)<0. (50) 



In equation l|47p . the projection P acts on the variables only; since / is a function of t,x,y,£, 
V(f) should be understood as 

P(f)(t,x,-)=P(f(t,x,-)), 
and the above equality holds between functions in L 2 (Y x R), almost everywhere in t, x. 

Proof. First step. Construction of f\. The existence and uniqueness of fx follows from a fixed point 
theorem in Xt- We define the application 4>t : Xy — > Xt by (frr(f) — 9, where g is the solution of the 
linear equation 

d t g + a(y,0-V x g + Xg = XP(f), 
g(t = 0, x, y, £) = u (x, y)) 

The existence and uniqueness of g follows from well-known results on the theory of linear transport 
equations (recall that a G C 1 ). Moreover, if /i, /a € and #j = (jyr(fi), i = 1,2, then g = gx — g 2 is a 
solution of 

d t g + a(y, ■ V x g + Xg = A [P{h) - P(f 2 )} , 
g(t = 0,x,y,0 = 0. 

Multiplying the above equation by g, and integrating on M.^ x Y x Kg, we obtain the estimate 

ollllflWII/^rR'VxV-xR) + A lbWlll^R-v x y xR ) < A / - "P(/ 2 )] g. 

Recall that the projection P is Lipschitz continuous with Lipschitz constant 1. Thus 



/ 



[P(h) - V(h)]g < \\\P(h{t)) - P(f2(t))\\l^xYxR) + JlbOHWxyxi 



xFx» 2 



< 2~ll(/l ~ •^ 2 )( i )Hi 2 (R N xYxR) + j HsWIll^R^xYxI 



Eventually, we obtain 



-n\\9(t)\\ 2 Li(R" X Y xR) + ^II5MIIl 2 (R«xYxR) < A||(/i - /2)(i)|lL2(R« x y xR ) < A||/i - h\\x T - 

A straightforward application of Gronwall's lemma yields 



\\9i-9 2 \\x T < Vl-e- AT ||/i-/ 2 ||x T . 
Thus (j>T is a contractant application and has a unique fixed point in X T , which we call f\. 
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Second step. L 2 estimate. Multiplying (|47|) by fx and integrating on R x Y x R, we infer 
\ d f 

nTll/AWIlL2(R»xFxl) + M\fx(t)\\ 2 L 2( R N xY xR) ^ A / V{f\)f\. 

Notice that 6 /C; thus the Lipschitz continuity of V entails that for almost every t, x 

Wfx)(t,x)\\ E < \\fx(t,x)\\ E . 

Hence, using the Cauchy-Schwartz inequality, we deduce that t i— > | |/a (*) | x y xR) ls nonincreasing 

on [Q,T]. The equality 

\x(t u o(x, y))\ 2 dxdyd£ = / \x(£,u Q (x,y))\ dx dy d£ = \u (x,y)\ dx dy 

'xYxR VxYxI JR N xY 

then yields the desired result. 

Third step. Compact support in £. Let us prove now that /a(-,£) = if |£| > M: let (p G £>(R) be an 
arbitrary test function such that <^(£) = when |£| < M. Then V{f\)<p = since V(f\) € IC, and thus 
/ A </? is a solution of 

^ (/ A p) + a • V, (/ aV j) + A (f x <p) = 0, 

{fx<p) (t = 0,x,y,t) = 0. 
Hence (/ A y?)(i, x, y, £) = for almost every t,x,y,£, and / A (-,£) = if |£| > M. 
Fourth step. Sign property. We now prove the sign property, namely 

sgn(C)/A = |/a | < 1 a.e. 

This relies on the following fact: if <? 6 1C, then sgn(£)g(j/, £) € [0,1] for almost every y,£. Indeed, 
£,(., £) = if £ < -M, and thus if —M < £ < 0, 



ff(»,0 = - / "(i/>£'K'<o. 

J-M 



Hence g(y, •) is non positive and non increasing on (— oo, 0). Similarly, g(y, •) is non negative and non 
decreasing on (0, oo). And if £ < < £', then 

g(y, £') - g(y, = 1 - / ^(y 5 «>) du> < 1. 



Hence the sign property is true for functions in /C. 
Multiplying (|47|) by sgn(£), we are led to 

§- t MO fx) + a{y, £) • V, (sgn(£)/ A ) + A (sgn(£)/ A ) = XV(fx) G [0, A]. 

And at time t = 0, sgn(£)/ A (t = 0) = |x(£, uo)| G [0, 1]. Thus, using a maximum principle for this linear 
transport equation, we deduce that the sign property is satisfied for fx- 

Fifth step. Uniform continuity at time t = 0. Let S > be arbitrary, and let /q := fo * x 9 s , with 5 a 
standard mollifier. Then fo(x) £ JC for all x £ R^, and thus / A — /q is a solution of the equation 

j t (fx - fo) + <V, • V, (fx ~ /o 5 ) + A (fx — /o 5 ) - A (V(fx) - P(/ 5 )) - o(», £) • (/o ** V0 4 ) . 

Multiply the above equation by (/ A — /q ) and integrate on R N x Y x R. Using once more the Lipschitz 
continuity of the projection V, we obtain 



2^ I | /A ~ /o I L 2 (R«xYxR) - ll a llL°°(rx(-Af,Af))ll/A - /o I Il 2 (R« x YxR) I l/o I U 2 (R« x YxR) 1 1 V6> 



II/a 


- f s \\ 


II/a 


Jo 1 1 



±\\f -f*\\ < <L 

mW JX J ° Ml 2 (R n xYxR) — § - 
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As a consequence, we obtain the following estimate, which holds for all t > 0, A > and 5 > 



s Ct s 

||/a(*) - fa | \l 2 (R n xYx$L) — ~ + ~ ^0 I L 2 (R«xYxK) " 

Hence the uniform continuity property is true, with 

W ® ^ l>o + 2 ||/o - /o|L 2(R « xV - ><K) ) • 

Sixth step. Inequality for M\. Inequality l|50p is merely a particular case of the inequality 

{V(f)-f,V(f)-g) E <0 
which holds for all / 6 E, for all g € 1C. 

4.2 The hydro dynamic limit 

In this subsection , we prove the following result : 

Proposition 3. Let (/a)a>o be the family of solutions of the relaxation model (|47j) . and let f{t) = u) 
be the unique solution of the limit system I|43p - lj46"|) with initial data x(£,jUo(x,y)). Then as A — > oo, 

fx -»■ / in L 2 ((0, T)xR N xY x E). 



□ 



The above Proposition relies on an inequality of the type 

if \fx-f\ 2 <rx(t), 



with rx(t) — > as A — > oo. The calculations are very similar to those of the contraction principle in the 
previous section; the only difference lies in the fact that f\ and / are not solutions of the same equation. 

Before tackling the proof itself, let us derive a few properties on the weak limit of the sequence f\. 
Since the sequence f\ is bounded in Xt C L 2 ((0, T) x M. n x Y x K.), we can extract a subsequence, which 
we relabel /a, and find a function g 6 £ 2 ((0, T) x 1^ xFxl) such that fx weakly converges to g in L 2 . 
Moreover, the sequence V(fx) is bounded in L 2 ((0,T) x x Y x E), for all T > 0. Hence, extracting 
a further subsequence if necessary, we can find a function h S L 2 ((0, T) x M. N xFxR) such that P(f\) 
weakly converges towards h as A — > oo. Notice that the convex set K. is closed for the weak topology in 
L 2 . Consequently, h(t,x) G K for almost every t,x. At last, 

V(fx) -fx = 

where the O is meant in the sense of distributions. Hence, g = h, and in particular, we deduce that 
g{t,x) £ K for almost every (t, x). 

We are now ready to prove the contraction inequality; consider a mollifying sequence 9 s as in the 
previous section, and set f s = / *t, x & 6 , fx = fx *t,x S 5 ■ Then 

dtf + a{y,0-\l x f 5 =M S , 
dtfx +a(y,0-V x f S x =M(. 

Let us multiply the first equation by sgn(£) — 2/| , the second by 2(/^ — f s ), and add the two identities 
thus obtained; setting F^' S ' = sgn(^)/' 5 + |/f | 2 - 2f s f*' , we have 

dtF 5 / + a(y,0 ■ V X F S / = M 5 (sgn(0 - 2/f ) + 2Mi(f x ' - f). 
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We integrate over (0, t) x 1^ x Y x K and obtain 

/ F 5 / (t, x, y, dy = [ [ M s (sgn(£) - 2/f ) 

VxYxI JO JR N xYxM v y 

+2/'/ M'M'-f*) 

Jo Jl^xYxl 

+ / F 5 /{t = Q,x,y,£)dxdydti. 

Jr n x7xl 

We now pass to the limit as 5' — > 0, with all the other parameters fixed. Notice that 



lim 



M{(fi-f 6 ) = 



JR N x7xl 



Xa(/a-/ 5 ) 



= -A 



< 0, 



o «"xy> 



(/a-^(/ a )) 2 



M(P(/a)-/ S ) 



JR N xYxI 



since x) G K. for all t, x. The passage to the limit in F^' S (t = 0) does not rise any difficulty because 
of the strong continuity of the functions f\ at time t — 0. Hence, we retrieve 

/ {(l/ tf (*)l-l/ 5 (*)l 2 ) + !/'(*) -A(*)l 2 } 

Jr n xYxR 

< f I M s (sgn(0 - 2/ A ) 

JO VxYxI 

+ / = 0)| - \.f(t = 0)| 2 ) + \f(t = 0) - x(^o)| 2 } , 

./R^xYxR 

and thus, integrating once again this inegality for t £ [0, T], 

r i -i/*i 2 ) -m 2 } 

JO Jr n x7xl 

< [ T dt\ f [ M 5 (s) (sgn(0 - 2/ A (s)) ds 

Jo Uo Jr n xYxR 

+T f {(\f 5 (t = 0)1 - \f 6 (t = 0)| 2 ) + \f(t = 0) - x(tu )\ 2 } . 

Jr n xYxR 

We now pass to the limit as A — > oo, with 6 > fixed. Then 

II/a _ / <5 || 2 2((0,T)xR«xYxR ^ lb _ / 5 || 2 2((0,T)xR™xyxEJ 

and 



lim / dt\ [ [ M S (s) (sgn(0 - 2/ A (s)) 

A^ooJo UO Jr n xYxR 



dt 



JRNxYxI 



M\s) (sgn(0 - 2 ff (s)) ds 



< 0. 



Thus, we obtain, for all 5 > 



lis-/' 



5 1(2 



li 2 ((0,T)xR Jv xyxI 



< T 



{(l/ 4 (* = 0)| - |/ 4 (t = 0)| 2 ) + \f(t = 0) - x(^>o)| 2 } • 

/R"xYxt 

We have already proved in the previous section that the family f s (t = 0) strongly converges towards 
X(£,«o) as S vanishes, due to the continuity assumption at time t = 0. Hence, we obtain in the limit 

I Iff _ f\\ 2 L 2 ((0,T)xR N xYxR ^ 0, 

and consequently, g — f. Hence the result is proved. 
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5 The separate case : identification of the limit problem 

This section is devoted to the proof of Proposition [TJ Thus we focus on the limit system in the case 
where the flux A can be written as 

My, = a o(y)g(0, witn div y a = o. 

The interest of this case lies in the special structure of the limit system; indeed, we shall prove that 
the function u, which is the two-scale limit of the sequence u £ , is the solution of the scalar conservation 
law (fT9| . In view of Theorem [TJ we wish to emphasize that Proposition [TJ implies in particular that the 
entropy solution of lfT9| satisfies the constraint equation 

div y (a (y)g(u(t,x;y))) 

for almost every (>0,i£ M. N ; this fact is not completely obvious when g ^ Id. We will prove in the 
sequel that u(t, x) actually belongs to the constraint space Ko for a.e. t, x. 

Before tackling the proof of the theorem, let us mention that the limit problem lfl9|) is not the one 
which is expected from a vanishing viscosity approach. Precisely, for any given 8 > 0, let uf be the 
solution of 

d t u £ 5 + div x A ( j, ttf J - eSA x u £ s = 0, 

with the initial data u £ s (t = 0, x) = uq (x,x/e). Then for all e > 0, u e s — > u £ as 8 — > 0; moreover, the 
behavior of u £ s as e — > is known for each 8 > (see [H [6]). In the divergence-free case, for all 8 > 0, 

]hnu e s =u(t,x) inLl ocl 

e— »0 

where u is the entropy solution of 

d t u + div x ((a) g(u)) = 0, 

with initial data u(t = 0, x) = (uo(x, •)). Hence, it could be expected that the limits s — » and 6^0 
can be commuted, that is 



which would entail 



lim lim ui = lim lim u% , 



lim u £ 

e->0 



In general, this equality is false, even in a weak sense: a generic counter-example is the one of shear 
flows (see for instance the calculations in |9j). In that case, we have N — 2 and A(y,£) — (ai{y2)£,0), 
and the equation l|19p becomes 

d t u -I- a 1 (y 2 )d Xl u = 0, 

with the initial condition u(t = 0,x,y) — Uo(xi,X2,y2)- It is then easily checked that in general, the 
average of u over Y is not the solution of the transport equation 

d t u + (ai) d Xl u = 0. 

We now turn to the proof of Proposition [TJ In view of Theorem [TJ it is sufficient to prove that the 
entropy solution of l[T9|) belongs to Ko for a.e. t, x, or in other words, that Ko is invariant by the semi- 
group associated to equation l[T9|) . We prove this result in the slightly more general context of kinetic 
solutions. The core of the proof lies in the following 

Proposition 4. Let uq e L 1 (R JV , L°°(Y)) such that Uo(x, •) £ Ko for almost every x € M. N . 
Let v = v(t,x;y) £ C([0, oo); L 1 (M Ar x Yj) be the kinetic solution of 

d t v(t,x;y)+div x (a (y)g{v(t,x;y))) = 0, t > 0, x G M. N , y eY, 
v(t = 0,x;y) = u (x,y), 

i.e. f l {t,x,y,^) := x(£, v(t, x; y)) is a solution in the sense of distributions of 

dtf'+Mv) ■ Vs/V(0 = # e m, t > 0, x € R N , y £ Y, £ £ M, 
/ 1 (t = 0, : r,y,0 = x(^^o(a ;i y)), [ * L) 

and m is a non-negative measure on [0, oo) x M. N x Y x M. 
Then for a.e. t>0,xeR N , u(t,i)eX , 
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Proof. First, let us recall (see that for a11 T > °> 

f 1 = lim f x inC^T];!, 1 ^ xF xR)), 

A— *oo 

where fx — f\(t 7 x,y,{;) (A > 0) is the unique solution of the system 

d t .h + My) • Vx/a g'(0 + a/a = A x (f, u x ), 

u\(t,x,y) = J R f\{t,x,y,£) d£, (52) 
/a(* = 0) = x«,«o). 

Moreover, for every A > 0, u\ is the unique fixed point of the contractant application 

C((0,T);L 1 (R N x Y)) -» C((0, T); L 1 (M jv x F)) 

where «2 = / and / is the solution of 

d t f + a (y) • v./ g'(0 + Xf = A X (£, ui), . . 

/(*=o)= X (e.«o). ( j 

Thus, the whole point is to prove that the space 

{u G C([0,T];i 1 (M Ar x y));u(t,a;) G K a.e} 

is invariant by the application 0a- 

First, let us stress that for all u G L 1 (F), 

ueKo div y (a(y)x(£, u)) = in 2>'(Y x R). (54) 

Indeed, if u G K , then for all 8 > 0, set us = u * 6 s , with 8 s a standard mollifier. The function us is a 
solution of 

divy(a u,5) = rs, 

and the remainder rs vanishes strongly in i 1 (F) (see the calculations in the previous sections). Since 
the function us is smooth, if G G C 1 (M Ar ), we have 

div y (a G(us)) = G'(us)r s . 

Passing to the limit as 5 vanishes, we infer div y (aoG(u)) = for all G G C 1 (R iv ). At last, taking a 
sequence of smooth functions approaching x(£j m )i we deduce that div y (aox(£,,u)) — in V' pei (Y x R). 
Conversely, assume that div y (aox(£,,u)) — 0; then integrating this equation with respect to £ yields 
u G Ko. Hence l|54p is proved. 

Now, let mi G C{[0,T];L 1 {R N x Y)) such that Ui(t,x) G K a.e. Then div(a x(£, «i) = 0). Let / be 
the solution of (|53|) : since oo G Ko, the distribution div y (ao/) satisfies the transport equation 

d t (div(o /)) + <?'(0«o(y) • V, (div(o /)) + Adiv(a /) = 0, 

and div(ao/)(t = 0) = because Uo(x) G Ko a.e. Hence dw y (aof) = 0; integrating this equation with 
respect to £ gives u-i G Ko a.e. 

Consequently, u\(t, x; ■) G Ko a.e. Passing to the limit, we deduce that v(t, x; •) G Ko a.e. 

□ 

Let us now re- write equation f5Tj) : setting b(y) = ao(y) — ao(y), we have 

dtf + <*oG/)V a /y(0 = S e m - ftfo)V x /V(0 =: Mi. 

If u G L°°{R N ), then v G L°°([0, oo) x R N x Y), and it is easily checked that f 1 and Mi satisfy the 
compact support assumptions. According to the above Proposition, f 1 also satisfies (|43|) . and thanks to 
the structure of the right-hand side, the distribution Mi satisfies f46|) . Thus f 1 is the unique solution 
of the limit system, and Proposition Q] is proved. 
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6 Further remarks on the notion of limit system 



Here, we have gathered, by way of conclusion, a few remarks around the limit evolution system introduced 
in definition [TJ The main idea behind this section is that the limit system is not unique (although its so- 
lution always is) , and thus several other relevant equations can be written instead of (fTUj) . Unfortunately, 
there does not seem to be any rule which would allow to decide between two limit systems. 

Let us illustrate these words by a first series of examples : assume that the flux is divergence free, 
and let 

N 

K := {/ G Ll c (Y x K), £ d Vi ( ai f) = in V}. 
i=l 

We denote by P the projection on K in L\ oc {Y x R). Precisely, consider the dynamical system X(t,y;£) 
defined by 

f X(t,y;Z)=a(X(t,y,Z),t),t>0 
\ X(t = 0,y;O=y. 

Then for all £ G R, the Lebesgue measure on Y is invariant by the semi-group X(t; £) because of the 
hypothesis div y a(y, £) = 0. Hence by the ergodic theorem, for all / G L\ 0C (Y x R), there exists a function 
in L\ 0C {Y x R), denoted by P{f){y, f), such that 

P(f)(y,0= I™ ^ f f(X(t,y;0,0dt, 

T^oo 1 J 

and the limit holds a.e. in y, £ and in Y x (— R, R) for all R > 0. 

Set a := P(a). Then if / is a solution of the limit system, / also satisfies 

dtf + a(y,0-V x f = M 

and /, M satisfy (J9J) and (HI]) - lO- Indeed, 

M = M + [a(y,0-a(y,t)}-V x f 

and the term [a(y, £) — a(y, £)] • V x (/* x </j)(t, x,y,£) belongs to IK- 1 for all t, x. Of course, uniqueness holds 
for this limit system (the proof is exactly the same as the one in section [3]), and thus this constitutes as 
legitimate a limit system as the one in definition [TJ In fact, in the separate case, Proposition [TJ indicates 
that the above system seems to be the relevant one, rather than the one in definition [TJ Notice that the 
distribution M satisfies the additional property 

M*t, x <P(t,x) eK 1 - yt,x. 

Let us now go a little further: let 9 € C 1 (R) such that < 9 < 1, and let 

afl(y,0 = »(0o(»,0 + (i-»(0)a(y,0- 

Then / is a solution of 

d t f + a {v>Z)-V*f = Mo, 

for some distribution Mg satisfying (fTi)) . Thus this still constitutes a limit system which has the same 
structure as the one of definition [TJ Hence the limit system is highly non unique, and it must be seen 
as a way of identifying the two-scale limit of the sequence f £ , rather than as a kinetic formulation of a 
given conservation law, for instance. We wish to emphasize that if the flux A is not "separated", that is, 
if the hypotheses of Proposition [TJ are not satisfied, then in general, the function u such that / = lj< u is 
a solution of the limit system, is different from the solution v — v(t, x, y) of the scalar conservation law 

d t v + div x A(y 7 v) = 0, 

where the flux A is such that d^Ai(y,^) — a,i(y,£). Indeed, the function v above is not a solution of the 
cell problem in general, even if v(t — 0) is. In other words, the set K is not invariant by the evolution 
equation 

iv 

4=1 
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where to is a non-negative measure and g — l£<„. 

Let us now assume that the flux A is not divergence free. Then there are cases where yet another 
notion of limit problem can be given: assume that there exists real numbers p% < P2, and a family 
{v(-,p)} Pl <p<p 2 , which satisfies the following properties: 

1. The function {y,p) 1— » v(y,p) belongs to L°°(Y x [pi,P2]); 

2. For all p £ [pi,P2]j v(-,p) is an entropy solution of the cell problem; in other words, there exists a 
nonnegative measure m(y, £;f>) such that f(y,£;p) = ^-i<v(y,p) is a solution of 



in; 



3. For all p £ [pi,p 2 ], (v(-,p)) Y = 0; 

4. The distribution d p v is a nonnegative function in L l {Y x [pi,j?2]); this implies in particular that 
for all couples (p,p') £ [pi,P2] 2 such that p > p' , for almost every y 

v{y,p) > v(y,p'). 

Under these conditions, one can construct a kinetic formulation for equation |T|), based on the family 
v(x/e,p) of stationary solutions of |T|), rather than on the family of Kruzkov's inequalities. This kind of 
construction was achieved in [7] in a parabolic setting, following an idea developed by Emmanuel Audusse 
and Benoit Perthame in [5] ; these authors define a new notion of entropy solutions for a heterogeneous 
conservation law in dimension one, based on the comparison with a family of stationary solutions. Let 
us explain briefly how the kinetic formulation for entropy solutions of ([I]) is derived: let u e be an entropy 
solution of H]). Define the distribution m e £ £>'((0, 00) x R N x (^1,^2)) by 

--{m («* -(?')), + w, (* (!•«') - A - (?•(;■')))] } • <*> 

Then according to the comparison principle (which was known by Kruzkhov, see [16t I17j). m £ is a 
nonnegative measure on (0, 00) x x [pi,p2]- Now, set 



f{t,x,p) := l v{hp)<uHt}X) £ L°°([0,oo) x M JV x [ Pl ,p 2 ]). 

Thanks to the regularity assumptions on the family v(-,p), we can differentiate equality |55|) (which is 
meant in the sense of distributions) with respect to p, and we are led to 



This equation is in fact the appropriate kinetic formulation in the heterogeneous case; its main advantage 
on the equation ([4]) is the absence of the highly oscillating term 



£ 



a-N+i 1 C<u e 
Notice that for all p £ [pi,P2] 5 

div„ (^^-a(y,v(y,P))) = in V' peT (Y). (57) 

This equation is derived by differentiating equation 

div y A(y,v(y,p)) = 
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with respect to p. Thus, if we set 



Ov (y , p) 

a(y,p) ■= — g- — a(y,v(y,p)), 



the vector field a G L X (Y x [pi,p2]) is divergence-free, and the same kind of limit system as in the 
divergence free cas can be made. Of course, the interest of such a construction lies in the simplicity of 
the structure of the limit system in the divergence free case. 

Definition 3. Let / G L°°([0, oo), L 1 (R Ar x Y x R)), u G L 1 n L°°{R N x Y). We say that f is 
a generalized kinetic solution of the limit problem associated with the family v(-,p) if there exists a 
distribution M. G ([0, oo) x l w x Y x R) such that f and M. satisfy the following properties: 

1. Compact support in p: there exists {p'xiP'x) £ [PiiP2] 2 such that p\ < p[ < p' 2 < p%, and 

SuppAf c [0, oo) x R N xYx \p'i,p' 2 ]; 
f(t,x,y,p) = 1 if pi <p <p[, f(t,x,y,p) = if p 2 < p < p 2 . 

2. Microscopic equation for f : f is a solution in the sense of distributions onYx (pi,f>2) of 

div y (a(y,p)f(t, x, y,p)) = 0. (58) 

3. Evolution equation: the couple (f,A4) is a solution in the sense of distributions on [0, oo) x R N x 
Y x (pi,p 2 ) of 

f d t (v p (y,p)f) + a(y,p) -VJ = M, ^ 
1 /(* = 0,x,y,p) = l v ( y ,p) <ao (x,y) =■ fo{x,y,p); 

In other words, for any test function (j> G V per ([0, oo) x R N x Y x (pi,P2)), 

f(t,x,y,p)v p (y,p) {d t (f)(t,x,y,p) +a(y,v(y,p)) ■ W x <p(t,x,y,p)} dtdxdyd£ = 
= - (<t>, M) v v , - I lv(y,p)<u (x,y)Vp(y,p)4>(t = 0, x, y,p) dx dy d£. 

JR N xYxR 

4- Conditions on f: there exists a nonnegative measure v G M* ([0, oo) x M. N xFx!) such that 

d P f = -v, (60) 
0<f(t,x,y,0 < 1 a.e., (61) 

I ll/( s ) ~ M\lhm"xYx( p1 , P2 ) ds ^ °- ( 62 ) 

5. Condition on A4: for all (p G X>([0, oo) x R N such that <p < 0, the function M *t,x <P belongs to 
C([0, oo) x R N ,L 1 (Y x R)), and 

Iyxm ( M **,* <P) (*> x > ') ^ ^ °' ({ft) 
Vip G Lf 0C {Y x R), div y (aip) = 0, and d^ip > 0. y ' 

We now state without proof a result analogue to Theorems [H [2] : 

Proposition 5. Let A G W^ loc (Y x R). Assume that a G C^ er .(r x R) and that a G W 1 ' 1 ^ X (jPi,Ps))- 
Let u Q G rfR" x Y) n L^R^Cpe^F)) smc/i tta* uofc, •) is an entropy solution of the cell problem 
for almost every x G R N . Assume furthermore that there exists p\ < p' 2 in (pi,p2) 2 such that 

v(y,Pi) < u (x,y) < v(y,p' 2 ), 

and let 

fo\X,y,P) '■= lv(y,p)<u (x,y) 

Then the following results hold : 
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1. There exists a unique generalized kinetic solution f of the limit problem associated with the family 
(v(-,p)) Pl < p < P2 with initial data /q. Moreover, there exists a function u £ _L°°([0,oo) x 1^ x Y) 
such that 

f(t,X,y,P) = lv(y,p)<u(t,x,y) d.e. 

2. Let u e £ L°°([Q, oo) x M. N ) be the entropy solution of (JXJ) with initial data uo(x,x/e). Let 
f(t,x,y,p) = 1-v(y,p)<u(t,x,y) be the unique solution of the limit problem. Then for all regularization 
kernels (p s of the form 



x e 



with (p £ V(R N ), J(p=l,0<(p<l, we have, for all compact K C [0, oo) x 



lim lim 

5^0 e^O 



u £ (t, x) - u * x (p s (t, x, ^ 



(64) 



Hence a whole variety of limit systems can be given, depending on the choice of the family of solutions 
of the cell problem. However, it is not obvious that any given system is "better" than another one. But 
the important result, as far as homogenization is concerned, is that all systems have a unique solution. 
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